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Time Response...cescsssssssccsseseescrsseeesesssssseessereseeees EEIO2 


Definitions 


The zeros of transfer function are the values of the laplace transform variable, s, that 


cause the transfer function to become zero. 


The poles of transfer function are the values of the laplace transform variable, s, that 


cause the transfer function to become infinity. 


G(s) = s+5 


s*+5s+6 


- The output respons: 


function. 


5 
First order system G (s) = e Second order system G (s) = 


First order system 


Mechanical system 


A yt) 


=| 
a a 
DIU + keg 





yu 
mm + by(t)= ult) 
di 


The order of a system is defined as being the highest 
differential equation, or being the highest power of s in the denominator of th 





zero:-5 , poles: -2, -3. 


© of a system is the sum of two responses: Forces Response and 
the Natural Response. The forced response is also called the steady-state response 
particular solution. The natural response is also called the homogeneous si 


OY 


lution . 


power of derivative in the 
e transfer 
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5 
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Electrical system 
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Fall 2012 A o age 


Scanned by CamScanner 








Second order system ae 
Mechanical system Electrical system 
R L 

ae oT) 
u(t u(t) vi) 

W T 

mo y = utn eE y Rc MY y = ur) 

at dt dt dt 


- The time response of a control system consists of two parts: the transient response 
and the steady state response. The transient response is defined as the part of the time 
response which goes from the initial state to the final state. The steady-state response is 


defined as the behavior of the system as t approaches infinity after the transients have 
died out. Thus the system response y(t) may be written as: 


y(t) = yr(t) + yssCt) 


Where y,(t) denoted transient response and yss(t) denoted steady state respons« 


A first-order system only has s to the power one in the denominator, while a second- 
order system has the highest power of s in the denominator being two. 


Types of the input functions ( or test input signals) commonly are: 
Impulse function - step function - Ramp function - Sinusoidal function. 
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Px: Given the sy stem in figure below, write the outpu:, C(t), in general terms. Specify the 
forced and natural parts of the solution. 
y 5) 
RG)=F [+2] Cl) 
rad E 


Pole zero map. 





System response (step response) 














_ +2) A, B _2/5, 3/5 
cts) 545) 5 545 F s45 
where 
_ (s +2) _2 
c= TS re 
pe 8 t2) E 
S tug 3 








Forced fesponse Natural response 
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i j llowing conclusions: 
From the development summarized in last figure, We draw the fo g 


is, the pole at the 
1. A pole of the input function generates the form of the forced response (that pole a 
origin generated a step function at the output). 


2. A pole of the transfer function generates the form of the natural response (that is, 
the pole at -5 generated e”). 


st s 
a : , . ‘i e 
A pole on the real axis generates an exponential response of the form e“, where -a is th 


pole location on the real axis. Thus, the farther to the left a pole is On the negative real axis, the 
faster the exponential transient response will decay to Zero. 


4. The zeros and poles generate the amplitudes for both the forced and natural 
responses . 


Exercise: What are the orders of the systems described by the following transfer functions: 








a) G(s) = ———— 
) ms +bs+k 
" b) Gís)= 
RCs+1 
c) Gis)= 


LCs* +RCs+1 


First order system 


The general form of a transfer function for a first-order system is 


a 1 
G(s) = — = 
, sta  TS+1 
If the input is a unit step, where R(s) = I's. 





an ey 
$ Ss 1 he 
Initial slope = aa A 


gpp 







63% of final value 
ar ¢ = one time constant 


eht 
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Taking the i 


oy) Y 


C(s) =Ris)G(s; == 


s(s+4) 


Taking the inverse transform, the step response is given by 








c(t) = epf) + ealt) =1- en 








When 1 = 1/a, 


—-al = 
e |ie =E = 0.37 


CU Menite =1— | 172 = 1 -0.37 = 0.63 
Transient Response Specification 


Time constant (T.) associated with this system is t = 1/a . The time constant tells how 


quickly the system responds. So, systems that respond quickly have large values of a, and 
systems that respond slowly have small values of a. 


Rise time Tr , is defined as the time for the waveform to go from 10% to 90% of its final 
value. 


2 pti 20 
T, = — 


+ — — — 


a a a 
Settling timeTs i~ úctined as the time for the response to reach, and stay within, 2% of its final 
value. 
cm 4 a 
g NL 


Steady state errorfss) The steady-state error is the error after the transient response has 
decayed leaving only the continuous response. The error signal: 


e(t) = r(t) — y(t) = 1- (1 — e7%t) = gal 


As t approaches infinity , 7% approaches zero and the steady state error is : 
d ess = e(%) = lim[r(t) — y(t)] =0 


Ex Consider a simple first order system G(s) = K/(s+a). Assume the unit step response viven 
in figure below. Determine the value of a , K and first order characteristic. 
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From the response, we measure the time constant, that is, the time for the amplitude to reach 


63% of its final value. Since the final value is about 0.72. The time constant is evaluated whe: 
the curve reaches 0.63 x 0.72 = 0.45. l 


From the curve when 





e(t) 50.45 A = Tc =0.13 second. 
Tc = — = 0.13 > as E 
| | > 0.13 
Rise time 
mm 22 _ 2.2 
r= > = 77 0.286 s 
Settling time 
T £TE= : > 








ER an E K R 
(an EA ES 
The transfer function | a 
ats) = s+77 
ooann E ÓEÓA< TT A o a 
CA a “Fall 2 REA 
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x 
Find the capacitor voltage in the network shown in figure below if the switch closes at t = 0. 


Assume zero initial conditions. Also find the time constant, rise time, and settling time for the 
capacitor voltage. 


À A 
W y 
122 | 
ia i 


0.5 F 7 


1 
Vas) _ Cs 




















F 5 i 10 e 5 
Using voltage division. Van...” RETE Since Vils)=< . Vets) = 9-2) s ` s+2 
(R=33) 
l 22 22 4 4 a 
Therefore vc(t) = 5 - Se-'. Also. T =5 . T= g > A 5 a 
Program: 
len=[1 Zi; 
=tf£ (num, de 
step (5*G) 
Computer response: 
Trensfe f1nc= 1085 
s +62 
Step Fesponse 
| 
o 
o 
= 
3 
1 
25 3 
ower a= ———_e_e— A Page 
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. 4 stems using MATLAp 
ae the time responses and find Tr Ts for the following sy AB 
G1 j= G2(s) = —, 163) = =e 
s)= 
e > s+10 
(a) 
MATLAB codel tia ai cd ÓN 
Poo > - System G1 
a | ' A a | System G1 Settling Time (sec): 1.3 
numl=3; _ E ' Fi Fise Time (sec): 0.733 
denl=[1 3]; a $: | , Ce ee 
= 05 =, 
Gl1=tf (numl,denl) = a 
subplot(3,1,1) + } 
step(Gl) ! 
title('(a)') of! —L = Las : po 
grid 0 02 04 0.6 0.8 1 12 1.4 1.6 1.8 
num2=5; Tea 
den2=[1 10]; 
G2=tf (num2,den2) ig l AS. - 7 
subplot(3,1,2) j p 
step (G2) 
title('(b)') 2 E 
grid = 04 pr --- ' System G2 
num3=5; E |! P E ida l Settling Time (sec): 0.391 
den3=[2 4]; 02} y ida wei 
G3=tf (num3 ,den3) Mis Ly 1% y 
subplot(3, 1,3) 0 0.1 02 0.3 04 05 0.6 
step(G3) Time (sec) 
title('(c)') 
grid pe (c) 
2 ee A 
O ee 
3 > | System G3 j System G3 -- 
Z ái + Fise Time (sec): 1.1 + S09 Time (sec): 1.96 
LS 0.5 af > : = - ae 
7 | |" 
0 3 Ga dl 
0 0.5 1 i 25 a 
«y 1] 
Time (sec) 
A.J.A Fall 2012 EER. 
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ZATLAB code2 
numl=3; 

denl=[1 315 

Gl= tf(numl,den1) 
num2=5; 

den2=[1 10]; 

G2= tf(num2,den2) 
num3=5; 

den3=[2 4]; 

G3= tf (num3 ,den3 ) 
step(G1,G2,G3) 
legend('G1 ', 'G2 
grid 


Amplitude 
b 


'+'G3') 


Step Response 
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unit step) of first order systems. 
o help with the calculations, Some poinN 


Ex: 

The following figures present t 
Determine the transfer function 
were given explicitly on the graph. 


he step responses (to a 


s of the two systems. T 




















Step Response R 
2 A A — ar s— 
j GEA System: tff | 
18 we Tre (sec): 25.1 | 
i Amplitude: 1.92 | 
System: tff mol | 
Tare (sec): 6.07 
Amplitude: 1.41 
| a 
/ 
: | / 
g | F 
2 | 
£ 08 / 
6 fi 
4 
i 
j 
SN y 
| | 
| 
' 
' Tire (sect 
Step response ofa system 
Step Response 
x | a Oooo 
see 7 | 
q . ve | 
ý Pi System x 
Tire (sec): 0.5 
System x Amplitude: 0.367 
Time (sec): 0.2 
Amplitude: 0.253 
o 
3 p? 
| | 
| 
| 
£ C4 o6 03 E. Ok 
Time (sec) 


Step response of a system 


Ü ae roe 
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Second order system 


A general form of a second order system is: 


d'y dy(t) 
— +20, —— + wo y(t) = ko- u(t) 
a E 


Transfer function: 





s` +2€0,5 +0; | 
k : the gain of the system 
X : the damping ratio of the system 


@,, : the (undamped) natural frequency of the system. 


The block diagram and signal flow graph representations of a second order system are given in 
figure below 





















5, 268 l 
R(s) 1 @ o, K Cs) 
1 
assume k-1, for simplicity oer 
Y(s) _ o; 
IR(s) s- +2G@,s+ o; 
The characteristic equation for the system is 
5° +2č0,S +0? =0 
Therefore, its roots are given by 
= + y = 1 
S1, S2 =—GQ, 20, 
S12 = 70, tjwny1- P= -oa tjwa 
Oa : Magnituce of the real part of the pole and is the exponential damping frequency. 
wq : Imaginary part of the pole and is called damped frequency of oscillation. 
A.J. A Fall 2012 Page | 
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? =pendin. 


Step response for the general second-order system 


K Kos + K3 
(Me o TT 
~ s(s2 + 2205 +05) S Ze 


as ing by partial 
where it is assumed that ¢< 1 (the underdamped case). Expanding by P 
fractions, 


f 2 
1 (8 + Zn) + bes on poe 
OS Erra 


`ch i i udent, 
Taking the inverse Laplace transform, which is left as an exercise for the st 
produces 


c(t) = 1 — etet (so wn V1 — ¿+ sin wn Y 1 — a) 


$ 
VET 








1 
E = 1] — — e ih cos(wn y 1 — ¢71 — q) 
! vie 
A 








where $ = tan (¢/ 4/1 = 2). 


A plot of this response appears in Figure below for various values of Z. We now see the 
relationship between the value of Ç and the type of response obtained: The lower the value of 


nse. The natural frequency is a time-axis scale factor and does 
use Other than to scale it in time. 


Ç, the more oscillatory the respo 
not affect the nature of the resni 
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pe 


nding pe ins ' 
g On the value of (X) three distinct cases preseni as shown in table below 


Three distinct cases of second order systems 


— Free distinct cases of second order systems =  žć  ć  <ž 


Case Damping ratio Roots of characte. .“*ic System 
value equation specification 
>i Sı + s, and both are real. Overdamped 
¿=1 Sı = $, and both are real. l Critically damped 

3 4 1 Sı and s} are complex, one being Underdampec 


the conjugate of the other. 





We see that the various cases of second-order response are a function of Č, they are 
summarized in figure below. 





e Poles Step response 
At) 
i 
Undaniped! 
jo s-plane c(t) 
t 
O<íf<l | | 
É 
Unuerdumped 





cl) 
s-plane 
g = | rod 
LOr 3 


Critically damped 





jo eu 
-yt 0. A 1 ns | ere 
¢>! oo — 0 | ia 
PR, yo- t iverdamped 
€ OB a 
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Response 








System Pole-zero plot 
i sl 
a) RES] so | ow 
Ya +b 
Genera 
c{?) cin = 14 0.1710 e 
if 1.1710 144 
CAs) 
J 
á di a 0.5 
5 +95 +9 Jas 
Ovecdamped 
LITA 
an dr = |e cose «4 sin Fr) 
Ld =1-1 67" costÑr-19.47") 
Os) : 


l 
(o) Ris 3 y cts) 
+ 25+9 


Underdampel 


CAs) 


l 
+2 


Undaijad 


G(s) 


1 : 
Hi | | aa) 
le) a = E > E 
s+65+9 


Crincally damped 


A 





LA 





cd) = lee 





a ee O 
O AS 


a NCW? «ral AN 
A A ESEE E OA A ——————— A — I'I-ÁAAQ_— E E TON eE 
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0<p<t: $1, = -to t joy/l- E (=tu,<0)  underdamped 
f=1; 5,9=-w critically damped 
>l 5.5 =-lo, +0 2-1 overdamped 
=0: s s= joy widamped 
C<0: 51.9 = -tot joyi- E (=f, >0) negatively damped 






iyi 


- Underdamped 








Vr tne 


c(t) 





0 0.5 1 15 2 2.5 3 3s 4 
Step responses for second-order system damping cases 





— = El Sar PE) _ _ __—eb 
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rales id Sia 
QO z= > 
o >” > BS 
M8 O 


eg: . ee H n 

Definition of transient response specificatio A 
rì 

úl Allowable tolerance 





1. Delay time Td: The delay time is the time required for the response to reach half the fina! 
value the very first time. Ta= (10-43) [un offic 


2. Rise time Tr: The time required for the waveform to go from 0.1 of the final value to 0.9 of 
the final value 


ae aes 
T. = 2167406 Tys apars 


@ny 1-8 wd wd A wn ds 





| 








T= “0 atan *(/1-2/7)  n-cos"U7) 
r 


3. Peak time Tp: The peak time is the time required for the response to reach the 


first peak of 
the overshoot. 


Tl Tt 
T, = = — 


w0.y1-E We 
4. Maximum (percent) overshoot Mp: The maximum overshoot is the maximum peak value 


of the response curve measured from unit. If the final steady-state value of the response differs 
from unity, then it’s common to use the maximum percent overshoot. 





Mp% = 08% = e @"/VI-#) 100 


Mp% = OS% = c(tp)—c(co) 100 = Cmax—Cfinal .100 


c(œ) l Cfinal 


5. Settling time Ts: The time required for the transient's dam 


ped oscillations to reach and stay 
within +2% of the steady-state value. 


4 4 

T; = = 

¿0 Oa 
: oe , ee 3 
For varying within 5% ot the final value, the setting time is i, = — 
Un 





E A PP EII > EVA neen A wen 
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Step responses of second-order underdamped systems as poles move: 
a. with constant rea] part 

ginary part 

mping ratio 


b. with constant ima 
C. With Constant da 


<ts) 





ia) 
ls) 


Frequency the same 


2 1 
— A «-plano 
i ss < -0 €F 
E 
-- ———== Pule 
se motion 
— e. x 








ds — 
. PR ra aia  _ Q _=>=>z 2 A br 
AA E e erre ge 
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Nao 
DE 
is 
Damping Ratio and Damping Factor f 


on the step response y(t) of the prototype 


R 7 d w ” . . . 
The effects of the system parameters & an n he roots of the characteristic equation. 


second-order system can be studied by referring to t 


i i 'ollowing 
The corresponding time responses for the following figure are obtained by the following 
sequence of MATLAB functions. 


cle 
clf 
w=10; 
for zeta=[ 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2] 

num = [w*2 i 3 Step Responce for diffrent yar.< of zeta 
den= [1 2*zeta*w w%2]; q 
t=0:0.01:2; 
step(num,den, t) 





hold on; 

disp('zeta and poles') 

zeta 

roots([1 2*zeta*w w%2]) 2 
end = 
xlabel('Time(sees)') =, 


ylabel( ‘Amplitude 


title( ‘Step Response for different values of 
zeta' ) 


grid 
mitin) toned 
clc as Step Fissponse for diferent vales of w 
clf 
zeta=0.7 ; 12 
for w=[ 5 10 15 20 25 ] 
num = [w^2 ] ; 
den= [1 2*zeta*w w%2]; y 
t=0:0.01:2; : 
step(num,den, t) 
hold on; 
disp('zeta and poles’) 
w 
roots([1 2*zeta*w w%2]) de 
end 





xlabel('Time') or 
ylabel('Amplitude' ) 

title( ‘Step Response for different values of w' ) 
grid 


AA __ NS A TEA eres 


AAA AAA AS meza 
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istea having the following 


pplied to a unity feedback s 
15 


EX A step input is a 
transfer function: 





G(s) =- 
i s(s + 3) 
522 mu 
s(s+3)+15 s* +3s-15 


Fi k 
ind (a) the closed-loop transfer function, (b) 0, (c) č and (d) a, 


open-loop 


Solution: 
(a) Closed-loop transfer function = G(s) = 
1+ G(s) 
(b) Comparing the denominator of the closed-loop transfer function with that of a 


standard second order transfer function, we get 
a, =15 
w, = V15 =3.873 rad/s 


or 
= ().387 





(c) From the same comparison, we also get 


TE 8 
(Als 
y NU) + 


(4 = 0, J1- 8? = 3571 rad/s 


contro! system 


(d) 
=x A second order control system, having ¿ = 04 and œ = 5 rad 
zusiected to a step input. Determine (a) closed-loop transfer function, (b) t,, (c) 
order co 


Solution: 
OA 


R(s) 5° +2%0,5+0% 


and (e) Mp 
(a) The standara cluseu-rup transie: iunction for 2 seeond 
C(s) _ : 

















given by 
Substituting the values of € and w,, we get 
R(s) s? +45+25 
+3 -1 m 
(b) jia oc 2-7 ae 4335 
n fi =e Se 
x x 
LS = 0.686 s 
$ Pa AE avi -016 j 
d Ln 
dd Eu, (ANS) 
esM, = 100 ex (_ ] E 
(e) Dl p SE DRR p Í l- ë? f 0.16 
A. LA p o LN Fall 2012 i 











\ 
= 100 exp 
j yin 
A en A 


1 
|=25.4 
} 
Pa 


, 
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: ine the value of 
EX In the closed-loop system shown in Figure determine | ies t 
K and a such that the de cent peak overshoot due to the unit step input does no 
exceed 10%. Settling time should be equal to 1 s. 





Fig. - Closed-loop system 
Solution: The overall transfer function of the system is 


190. _£ 


R(s) s(st+a)+K 


K 


s +as+K 


Comparing the characteristic equation st 2 as + K = 0, with that of a standard second 
order system, we gel 








ay = JE 
ae s ti 
20% 2VK 
4 a 
Now, t= > =1s(given: 
50 
E, =4 


= A (using Eq. (Mm 
f 8 


or E 





E (10 
The peak overshoot condition gives 
f » 
exp | - xe | = 0.1 
L h-E) 
= — ates = 2.303 (iii 
ui go 


Substituting the valve of € from Eq. (i) in Ec. (iii), we get 


Ra 
ae 
Ra = 2.503 
a 
4 am oi 
y 4K 


ao m2303 
Jak - e 


Putting Une value of g from Eq. (ii) in this equation, we ge: 





28) __2303 
Jak —64 
a E ae 
2.303" 
> 1.6 sl ange 
o A a 1 ` 
J.A Fall 2012 O Page 
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EA For the second-0; 


der systems that follow, find E, a. TE 1P,15,0S. 





À 16 0.04 
e E ran 


(a) 


5 
On| = 16 r/s, 


20 = 3. Therefore 5 = 0375, 











On =4 
4 fe ares 2 Step Pesporse 
— Su ys Overshoot (%) 28 1 
Ts = Con sa 2.667 S At ire (sec) 0.85 
x 6 
Tp= = 0.8472 5 << 
rd aning Tema (sec) 2 66 
On 1-¢ NG e 
> af a l 
OS =e o V1-C 100=28 sung 
has e (sec) 0356 
. 2.167+06 2.16+0.37540.6 F 
T; = — - = — ir 
0, 4 
T, == 0.35 a 
CLE 
num=16; 
Ł=0:0 91:10, = 
den=(1 3 16]; 


G=tf(num,den) ; 

omegan=sqrt(den(3)); 

zeta=den(2)/(2*omegan); 

Ts=4/(zeta*omegan); 

Tp=pi / (omegan*sqrt (1-zeta"2)); 
OS=exp(-zeta*pi/sqrt(1l-zeta”2))*100; 
Tra=(1.76*zetaa”3 - 0.417*zetaa*2 + 1.039*zctaa 


+ 1)/omegana 
[y,t] = step(G,t); 


step(G,t) 
disp(' «mega zeta Tr : Ts OS%') 
{ omegan zet: TE Tp Ts OS] 
ome ge zeta Tr Tp Ts 0S% 
4.0000 0.3750 C.3600 0.8472 2.6667 0597 
SL I A A A A eT a ee a Sere 
LA Fall 2012 Pag: 
j 
j 
= ATES BEE Er hon ne the rtf 





a Ge 
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2 another solution 

cic 

num=16; 

t=0:0.01:5; 

den=[1 3 16]; 

G=tf (num, den); 
omegan=sqrt (den(3)); 
zeta=den(2)/(2*omegan) ; 
[y,t] = step(G,t); 

yl= find(y>=0.1); 

y9= find(y>=0.9); 

ys= find(y>0.98 £ y<1.02); 
yd= find(y>=0.5); 

n= max(y) 

yp= find(y>=n) ; 

Td= t(yd(1)) 
Tr=t(y9(1))-t(y1(1)) 
Tp=t (yp(1)) 

Ts=t (ys(45)) 


oP 90 90 ye 





time delay 
rise time 
peak time 
settling time 


os = (max(y) -y(end)/y(end))*100 ¿ overshoot 
step(G,t) 
disp(' omega zeta Td Tr Tp Ts OS%') 
[ omegan zeta Td Tr Tp» TS OS] 
System G 
Feak amplitude: 1.28 | 
14 Overshoot (%): 28.1 Sree 
At time (sec): 0.25 
1.2 7? Y 
/ ' System G 
| Settling Time (sec): 2.66 
¿SSPMRGz_— Zo . P Pope 
o J Fase Time (secj0386 -~ OT 
T 
2 ve rl | 
2 
a 
Rutt 1 
3.5 
Time (se: 
Tr Tr Ts 05% 
36 0.85 2.671 28.059 
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System, Response with additional Poles 


If a system has e 
he es as more than two poles or has ZETOS, We cannot use the formulae to calculate the 
p ance specifications that we derived. 


Under certain conditions, a s 


: ystem with more than two poles or with zeros can be 
approximated as a second-or 


der system that has just two complex dominant poles. 


Once we justify this approximation, the formulas f 


. i or percent overshoot, settling 
time, and peak time can be applied to these higher 


-Order systems. 


Consider a three-pole system with complex poles and a third pole on the real axis. Assuming 


that the complex poles are at 70, +j0n/1— 72 and the real pole is at a,. The step 
response of the system in time domain. 


c(t) = Au(t) + e (B cos wat + C sin wat) + De 


If |a,| >> |w,,| the pure exponential will die out much more rapidly than the second-order 
underdamped step response. If the pure exponential term decays to an insignificant value at 
the time of the first overshoot, such parameters as percent overshoot, settling time, and 
peak time will be generated by the second-order underdamped step response component. 
Thus, the total response will approach that of a pure second-order sysiem. 


We can show, through a partial fraction expansion, that the residue of the third pole, ina 
three-pole system with dominant second-order poles and no zeros, will actually decrease in 


magnitude as iis third pole is moved fariher into the left half-plane 
Ex 


num=5; 
den=poly([ -1 -2 -20]) 
[r,p,k]= residue(num,den) 


r= p= 
0.0146 -20 
-0.2778 E 
0.2632 


-1 


- If the real pole is five times farther to the left than the dominant poles, |a,| > 5!gw 


nl 
we assume that the system is represented by its dominant second-order pair of poles. The 
step response of the system in time domain. 


c(t) = Ault) +e “7 {B cos wat + Csin wat) 


_ _—.— A > 6 ER 9 AVE REE ET AIO AAA AAA AA EEE i 
A.J. A 
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EX 





Determine the vaii lity of « second-order appr m for each of these two transfer functions: 
700 36 
a. G(s) ZAR b. GL 350 
(s+ 15)(s + 45 + 100) (5+ 4)? +25 + 90) 


The second-order approximation is valid, since the dominant poles have a real part of 
-2 and the higher-order pole is at —15, i.e. more than five-times further. 
b. The second-order approximation is not valid. since the dominant poles have « real part 


of —1 and the higher-order pole is at A. i.e. not more than five-times further. 


numl= 700; 

denl=conv([{1 15],[1 4 100]); 
polesl=roots(den1l) 

gl= tf(numl,denl)>; 


num2= 360; eae 
denz=conv(,i 4],[1 Z 90]); i iig = 
poles2=roots(den2) L—e 
g2= :f(num2,den2); 


step(gl,g2) ESTAFA 


Amplitude 
€ 


poles1 = 


-15.0000 za% 
-2.0000 + 9.7980 Ai y 
-2.0000 - 9.7980; | 


| 
poles2 = J 


-1 .0000 + 9.4340i Tere (sec) 
-1.0000 - 9.4340i 


-4.0000 











E 
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AX Find the Step response of 
Se of each e follow Sfer f i i 
at of the following transfer functions using 


b- Simulink c- GUI 
24,5 2 
Wap TM ne OH 
+ 24.542 GHIO + 454+ 24542)" GEI 445 $ 24342 ) 
a- MATLAB l 


numl =[24.542]; 

denl=[ 1 4 24.542] 

Tl= tf (numl,den1) + aap 
num2 =[245.42]; . ~~ a 
den2=conv({1 10},(1 4 24.542)) 1 ’ S 
T2= tf (num2,den2) 

num3 =[73.626]; 

den3=conv([1 3],[1 4 24.542]) 
T3= tf (num3,den3) 
step(T1,T2,T3) 
Llagend( "11", TZ", 23") 


Step Sesponse 


$ 
E 
2 


b- simulink 





24.542 
+45+24.542 






246.42 





A e > — 
conví[1 10].(1 424.542)X5) 4 


Scope 








73.626 


confi 3).[1 4 24.542]Xs) 
13 





c- GUI : Another method to obtain systems response is though the use ltiview 
>> Itiview 





1- Select file 2- select import from file menu 
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System Response with Zeros 


ritten 1 > of the two general forms 
The transfer function for a second-order system Can be written in one g 


(depending on whether the system has a Zero OF not). 


K(s+a) _ G1(s 
Case1 : Gifts) = a eee `. Case2: G2(s) = lose = sG1(s) + aci ) 


s? +bs +c 


; , 29 e AO -3, -5, and -10. 
Starting with a two-pole system with poles at (-1 +j2.828), we conseculs ely add zeros at -5, -3, 


a 9 3(s + 3) 1.8(s + 5) e 0.9(s + 10) 
P= EA ta 4 Os+9' METAS 1° ste 2849 


The results, normalized to the steady-state value, are plotted in Figure below. 


step Response 
deng=[1 2 9]; 

Ta=tf ([1 3]*9/3,deng); 
Tb=tf([1 5] *9/5,deng); 
Tc=tf ([1 10] * 9/10, deng); 


T=tf (9, deng); A A naam namie’ 
step(T,Ta,Tb, Tc) fw floes 
text (0.5, 0.6, ‘no zero') E lÍ zero at -10 
text(0.4,0.7,'zero at -10') oe [no zero 


text(0.35, 0.8, ‘zero at -5') 
text (0.3,0.9,'zero at -3') 


Tare ¡sec 


- The motion of Case 2 systems is complicated somewhat over that of Case 1 systems 
by the presence of the zero ( s+a) in the transfer function. 


- The zero has little effect on the settling time of the system, but can significantly affect the 
overshoot. 


- How much effect the zero has on overshoot depends on where it is (along the real 
axis) compared to the poles of the system. 


o Zero is far to the left of the poles: it has little affect on the motion of the system. This 
makes the motion much like that of a Case | system. 


o Zero located near or inside the poles: it wi!) significantly attect the overshoot of the system 








AE N SAREE A E E SEU A) 


AJA Fal! 2012 Page 151 


Scanned by CamScanner 






t valuating Pole Z y \ e € 
` 7 ery Cancellati ` 

« ion l SINO \ Sid 
as 4 u > 


TO as shown 
T(s) = k(s + 2) 
(S+p3)(s? + as +b) 


If the > te 
he pole te 3) , and the Zero term(s+z 


rm, (s + p 


), Cancel out. 


We are left with 


T(s) —__Kis+z) 

3) (5 + as +5) 
aS a Second-order transfer function. From 
the pole at =P3 , the 


another perspective. ji 
n a partial-fraction e 
exponential dec 











the zero at - z is y ery close t 
Xpansion of T(s) wil] Show that the residue f the 
ay ts much smaller than the amplitude of the second-order response. Let us 
look at an example. 
Ex 
Determine the Validity of a second-order Step-response approxima-tion for ich transfer 
function shown below, find the approximate response. 
26 26.25(s +4) 
Ci Ete) O Cia) = ee = = 
S($+3.S)(s+5)(s— 6) Sis + 4.01)(s + S)is — 6) 
The partial-fraction expansion of Cl(s) 
Lo 3 33 i 
lee are i s E s4a8 
The residue of the pole at -3.5. which is Closest to the zero at -4. is equal to? and is not 
negligible compared to the other residues, Thus, a second-order Step response 
approximation cannot be made for C1(s). 
The partial-fraction expansion of C2(s) 
0.87 53 44 0.033 
TE ares ares 
The residue of the pole at -4.01, which is closest to the zero at -4, is equal to 0.03 4 about ty 
orders of magnitude below any of the other residues. Hence, we oo a nd-ordej 
approximation by ne glecting the response generated by the pole at -4.01: 
, 0.87 53 44 
ey Se San 
and the response c>(1) is approximately 
C2(t) = 0.87 — $.3e-* + 4.400 
EV O Skall 2012 
A.J.A 


TIS. A ananass o anne 
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Laplace Transform Solution of State Equations 


Consider the state equation 
: x = Ax + Bu 
and the output equation 
y = Cx + Du o 
Taking the Laplace transform of both sides of the state equation yields 
Sx) - 70) = AXG) + BUCS) 
sX(s) -x(0) = A(s) +BU(s) 
A 


(sI -A) X(s) = x(0) +BU(s) 


X(s) = (sI -A)” [x(0) +BU(s)] 
Taking the Laplace transform of the output equation yields 
Y (s) = CX(s) + DU(s) 
Note : The eigenvalues of the system matrix A = Roots of det(sI-A)=0. 
Ex: 
Given the system represented in state space 


2 
=| 0 el, pe x0)=| | y=[1 3]x 
3 -5 1j 1 


do the following: 


a. Solve for y(t) using state-space and Laplace transform techniques. 
Ł Tind the eigenvalues and the system peles. 











s —2 | ] s+5 2 
a. Since sI- A = ll (sl-AY == . Also. 
3 s+5| s+55+6| -3 s 
0 
BU(s)= š 
1/(s+1) 
4 | (As +7547) 
The state vector is X(s) = (s1 -AY [x(0)+ BU(s)] =——- —— 5 i 
(s+lXs+2is+3)| s—-4s-6 
Ss- +25-4 i5 2 5 
The output is Y(s)=[1 3]X(s) = Me =- a i2 RS 
(s+1)(s+2)(s+3) st] s+2 543 
Taking the inverse Laplace transform yields v(1)=-0.5e * -12e + CESA”. 


b. The eigenvalues are given by the roots of [sI — Al = s` +55+6. or 2 and -3. 
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MATLAB code 
cle 
syms s 
A=[0 2;-3 -5]; 
B= 0517; 35 exp(-3 t) pb 
Ee Be eS > 2 Jexpi=c 0) = AAA 
XO=[2;1]; i 
U=1/(s+1) 
I=eye(2); 
X= ((S*I-A)*-1)* (X0+B*U); 
x1l=ilaplace(x(1)); 
x2=ilaplace(x(2)); 
Y=x1+3*x2; 


> “u(t)=exp(-t) 


Eigen valus 


E= 
pretty(y) 
'Eigen valus' -2 . 0000 


E=eig(A) -3 . 0000 


Time Domain Solution of State Equations 
The solution in the time domain is given directly by 


EA E 
x(1) = eNx(0) + | ex Bu(rjdz Tp): 1 


0 


| = D(1)x(0) + [ou — )Bu(r)dr | 


Q(t): state transition matrix, 0(t) = e^! 


wns o) = A 


ES —A) | 


| Dl) = 2a - A)! = "| e 
| (sI — A) | 
EX 


For the following state equation and initial state vector. Find the state-transition matrix and 
then solve for x(t). Where u(t) is a unit step. 


a E TE... 10 
“=| 0+ y 


PA A A o 0 ETS A A E A 
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We use the fact that Qíi) is the inverse Laplace transform 
‘sl — AY. Thus, first find (s1 — A) as 


I 
a-a- (s +6) 





from which 
F -8 5 +6548 s*+65+8 
lI- a = 
ita 2 +4+654+8 -8 sS 


?+65+8 s.4+65+8 


Expanding each term in the matrix on the right by partial fractions vields 
2 1 1/2 _ 1/2 ) 
s+2 5+4 s4+2 s+4 


—4 4 —] Z 2 
43 sta s+2 s+4 


Finally, taking the inverse Laplace transform of each term, we obtain 


(2e-* a e") Ga e ze“) 


(sI - A)” = 


Db(1) = 2 2 
(—4e-* + 407) (—e7™ + 2e) 
1 ~2(1-1)} l Al! ») 
=e -xe 
P(r- 1)B = | ( 2 | 
i(- ein + 2e ute) ] 
AU _ A 
@(1)x(0) = | pener | 
(se? + 4e*) 


La [ea -les T e“dr 
2 o 2 0 


t t 
-e> f erdre | edr 
0 0 


ll 


/ “(1 - 1)Bu(ejdr 
0 


1_1 -u plot 


8 4 8 
ly 1 
- SAS y 
x(t) = (x10) + | M(t - 1)Bu(r)jdr = 
0 1 i 7-4 
2 A — AAA A A A A 
A.J.A Fall 2012 Page 1: 
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st ATLAB code 

cic 

syms s t tau 
A=[0 1;-8 -6]; 
B=[ 0;1]; 
xXO=[1;0]; 

U=1 ; Su(t)=step function 
I=eye(2); 

E=((s*I-A)”-1); 

Fi= ilaplace(E); 


pretty(Fi) 
Fitmtau=subs(Fi,t,t-tau); 
"Pi(t-ta)" 
pretty(Fitmtau) 
x= Fi*XO+int(Fitmtau*B*1,tau,0,t); 
' x (t) ' 
pretty (x) 
+- a + 
| exp(-2 t exp(-4 t) | 
| p(-2 t) - exp(-4 Tt), -=- | 
Fi= | 2 | 
| | 
| 4 exp(-4 t) 4 exp(-2 t), 2 exp(-4 t) exp(-2 t) | 
+- - 
-+ 
exp(2 tau - 2 t) exp(4 tau - 4 t) | 
| 2 exp(2 tau - 2 t) - apta tai 4 Uh == eS=ecotón pr essen | 
> ? 
Fitt-ta) ~ | 2 2 
| 
| 4 exp(4 tau - 4 t) - 4 exp(2 tau 2 t), 2 exp(4 tau - 4 t) - exp(2 tau - 2 t) | 
+ + 
+- -+4 
| 7 exp(-2 t) 7 exp(-4 t) | 
[ ess > E SNS + 1/8 | 
| 4 8 | 
x(t) = | | 
| 7 exp(-4 t) ? exp(-2 t} | 
| po. SSS E O | 
| 2 | 
+- + 
E ES Re SUS A E wee ADA SRE 1 AND AG a smi A SS i nur mar — 
AJA Fall 2012 Page 15 
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B dl > 
` 


x ynse a roxima-tion for each transfer { action 
] ralidit econd-order step-response ppro- 
Determine the val id MO of a secon 


shown below. 


197.14(s +7) 
185.71(s +7) b. G(s) = MT 
a. G(s) = Gy G5)(s + 10)(s + 20) ©) = +69) + 10) +20) 


1 0.8942 15918 0.3023 
a. Expanding G(s) by partial fractions yields G(s) =- + 5320 TT PFE 

A 
But -0.3023 is not an order of magnitude less than residues of second-order terms (term 2 


and 3). Therefore. a second-order approximation is not valid. 
0.9782 1.90/70 9078 0.0704 0704 
54+20 5410. $465 





l 
b. Expanding G(s) by partial fractions yields G(s) == + 
5 


. . : > > 2 
But 0.0704 is an order of magnitude less than residues of second-order terms (term - and 


3). Therefore, a second-order approximation is valid. 





Ex 
For the system shown in Figure below, a step torque is applied at 6,(1). + ‘nd 
a. The transfer function, G(s) = @2(t) / T(t) b. The percent overshoot, settling time, and peak time. 
TN 8n iry 
A A \ 
C+ E 
1 N-m-s/rad 
Writing the loop equations, (s*+s)a,@)-s 0,()=7 (5) 


-s 0, (s)+(s -1)0,(5)=0 


pas T(s) 
- -3 0 2 T (s 


Solving for 02(s). 





ri (s) 
sts y S “+s +] 
-s 1 

Forming the transfer function. 

] 

ols) 

T (s) si+s+1 
Ihus op = 1,230) = 1. rhus, ¢ = 0.5. %OS = 16.3%. Ts =8 seconds Tp = 3.63 seconds. 


-vwy a ee o E: waes 
J.A 
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tx: For the s- 


velocity Of th 


«a 
~ 


a % è 
TM Shown in fi 


| figure E 
mass to M: f Eure below 









“l; (5) find an equation that relates nio. oo Telates settling time of the 
= ates rise time of the Velocity of the mass to M. 
Y nung the equation of Motion 
i (Ms +8 =} 
Thus. the transfer function is. ia iad, “i 
X(s) | 
yo neee 
F(s) Mi + 8s 
Difterentiati yield 
tte ‘ME to yield the transfer function in terms of velocity. 
sX(s) l I/M 
T es = Ā~— 
Fis) Ms 4g 142, 
M 
Thus. the seni 


ing time, T,. and the rise time. 7, are given by 
4 
T, -> 

€ KIM 


N AS 
M: === = .275) 
"SM 0.275M 


t3 | 


Ex For each of the transfer functions sho 
plot them on the S.-plane, and then write 
response without solving for the inv 
(overdamped, underdamr 


wn below, find the locations of the poles and zeros, 
an expression for the general form of the step 

erse Laplace transform. State the nature of each response 
ed, and so on). 


de 5 
a. T(s) = 





— 3 ah 
— 5 — 
ae ia (s+ 3)(s3 6) 
10(s + 7) 
te Caza San 
20 
d. Tís) 


~ s? + 6s — 144 
pi AR 
e. si! dak ar 
pm o A 5) 
> PAS (s + 10)* 
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a. Pole: -2; e(t) = A + Ber?t; first-order response. 


b. Poles: -3. -6: c(t) = A + Be! + Cet; overdamped response. 





e. Poles: -10, -20; Zero: -7; c(t) = A + Belt + Ce-20t, overdamped response. 
d. Poles: (-3+j3/13 ), (3-j3/13 ) ; e(t) = A + Be3t cos (34/15 t + 4): underdamped. 
e. Poles: j3, -j3; Zero: -2. c(t)= A +B cos (3t + q); undamped 


f. Poles: -10, -10: Zero: -5; c(t) = A + Be! + Cte! critically damped. 


Ex Write the general form of the capacitor voltage for the electr cal network shown in Figure 
below. 


R, = 10kQ 


vtr) = utr} 





Writing the node equation at the capacitor, 


Vols) - V(s 

Vets) RE +E r Cs) A = 

EN 

Vels) R¡ 
tence, gy ee a e 
Ri + R> + Ls* Cs 

10 

s*+20s+S00 


as 10s 
52+205+500 


The step response is 


The poles are at 
-10 +20. Therefore. velt) = Ae!" cos (20t + ò). 
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Problem 1 HW 5 


A RC circuit 5 
has the following transfer function: Y(S) _ 2 For a step i > 
: =——— Fora ste =2V 
what is the time taken for th R(s) 10s+4 Ep InP RG ac 
© output of the RC circuit to reach 95% of its steady-state response? 


Problem 2 


A system has transfer function: 14S) _ 50 


R(s) s+50 


Find the time ; Wo da 
time constant. T. the Settling time, T.. and the rise time, T, for a unit-step input. 


Problem 3 


Consider the unit-step response of a unity-feedb 


5 noun ack control system whose open-loop transfer 
function is 


> l 
E ST 


Obtain the rise time, peak time, maximum overshoot, and settling time. 


Problem 4 
A second-order system is underdamped with a damping ratio of 01 and a natural frequency of 10Hz 
Find: 


a) the transfer function 
b) the time response when it is subjected to a unit-step input 


d) the rise time 


Problem 5 


Given the transfer function: 
100 
G(s) = ——______ 
s +15s +100 


find T,. PO, T, and T, 


Problem 6 
Obtain the unit-impulse response and the unit-step response 0! a unity feedback system whos: 
open-loop transfer function is 


G(s} = ZS 
s 
EI ES IIS 2 ES ET AA A A A A et e r E- 
A.J. A Fall 2012 Page ! 
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Problem 7 
The following block diagram of a space-vehicle attitude control system. Assume the tim 
constant T to be 3 sec, and the ratio K/J to be (2/9) rad’/sec. Find the daming ratio of the 


system. 












Space 
vehicle 


Problem 8 


Obtain the unit-step response, unit-ramp response, and unit-impulse response of the following 
System and repeat your solution using MATLAB. 


C(s) 10 





RS s$=25+10 


Problem 9 


Obtain the unit-step response, unit-ramp response, and unit- -impulse response of the following 
system and repeat your solution using MATLAB. 


HERRING. 
of * | 


where u is the input and y is the output. 


Problem 10 


Given the following state-space representation of a system, find Y(s): 


i E al $ [7 frina: 


y=[1 3]x; x(0) = g 





=  — ——_—_ _ __—_=3x-  __— e —— 4 € a 
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EXPERIMENT [BS A eaccoononnonooonocninnooonanionronocioronroo EOL 
OBJECTIVE 


l- Evaluate the effect of pole and zero location upon the time response of first- and 
second-order systems. 


2- Evaluate the effect of additional poles and zeros upon the time response of 
second-order systems. 


MINIMUM REQUIRED SOFTWARE PACKAGES MATLAB, Simulink, and 
the Control System Toolbox. 


PART 1 
PRELAB 





, ee. a E ae 
l- Given the transfer function G(s) = e evaluate settling time and rise time for 
s+a 


the following values of a: 1, 2, 3, 4. Also, plot the poles. 





2- Given the transfer function G(s) = =——~. 
si+as+b 


a- Evaluate percent overshoot, settling time, peak time, and rise time for the 
following values: a = 4, b = 25. Also, plot the poles. 


b- Calculate the vaiues ot a and b so that the imaginary part of the poles remains 
the same but the real part is increased two times over that of Prelab 2a, and 
repeat Prelab 2a. 


c. Calculate the values of a and b so that the imaginary part of the poles remains 


the same but the real part is decreased (1/2) time over that of Prelab 2a, and repeat 
Prelab 2a. 


3. a. For the system of Prelab 2a, calculate the values of a and b so that the real part 


of the poles remains the same but the imaginary part is increased two times over that 
of Prelab 2a, and repeat Prelab 2a. 


b. For the system of Preiab 2a, calculate the values of a and b so that the real part 
of the poles remains the same but the imaginary part is increased four times 


over that of Prelab 2a, and repeat Prelab 2a. 


TAr tha e ñ p “i a | E 1 aint y ' . 
4. a. For the syuicin ot Prelab 2- . calculate th: vaiuics of a and b su that the damping 
pS 
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j ency is increased two times over 
ratio remains the same but the natura! frequency 1s Increas 
that of Prelab 2a, and repeat Prelab 2a. 





b. For the system of Prelab 2a, calculate the values of a and i bos od ee 
ratio remains the same but the natural frequency is increased four ti 
that of Prelab 2a, and repeat Prelab 2a. 


5. Briefly describe the effects on the time response as the poles are changed in each 
of Prelab 2, 3, and 4. 


LAB 


|. Using Simulink, set up the systems of Prelab 1 and plot the step response of 
each of the four transfer functions on a single graph by using the Simulink LTI 


Viewer. Also, record the values of settling time and rise time for each step 
response. 


2. Using Simulink, set up the systems of Prelab 2. 
plot the step response of each of the three transfer 


Also, record the values of percent overshoot, settli 
time for each step response. 


Using the Simulink LTI Viewer. 
functions on a single graph. 
ng time, peak time, and rise 


3. Using Simulink, set up the systems of Prelab 2a 
LTI Viewer, plot the step response of each of the th 
single graph. Also, record the values of percent ove 
time, and rise time for each step response. 


and Prelab 3. Using the Simulink 
ree transfer functions on a 
rshoot, settling time, peak 


POSTLAB 


1. For the first-order systems, make a table of calculated and experimental values of 
settling time, rise time, and pole location. 

2. For the second-order systems of Prelab 2 
experimental values of percent overshoot, 
location. 


, make a table of calculated and 


settling time, peak time, rise time, and pole 


3. For the second-order systems of 
calculated and experimental values 
time, rise time, and pole location. 


Prelab 2a and Prelab 3, make a table of 
of percent overshoot, settling time, peak 
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bs nant har pi systenis of Prelab 2a and Prelab 4, make a table of 
calculated and experimental values of percent overshoot, settling time, peak 


time, rise time, and pole location. 


5. Discuss the effects of pole location upon the time response for both first- and 


second-order Systems. Discuss any discrepancies between your calculated and 
experimental values. 


PART 2 
PRELAB 


1. a. Given the transfer function G(s) = evaluate the percent over- 


5 
l s?+4s+25 
shoot, settling time, peak time, and rise time. Also, plot the poles. 


b. Add a pole at -200 to the system of Prelab la. Estimate whether the transient 
response in Prelab la will be appreciably affected. 


c. Repeat Prelab lb with the pole successively placed at -20, - 10, and -2 


2. A zero is added to the system of Prelab la at -200 and then moved to 5), - 20, 
- 10, - 5, and -2. Lis: the values of zero location in the order of the greatest to the 
least effect upon the pure second-order transient response. 


(Gera let a=3 and b= 3.01, 3.1 
A a 


3.3, 3.5, and 4.0. Which values of b will have minimal effec: upon the pure second- 
order transient response? 


3- Given the transfer function G(s) = 





Ejea l 30 and b 
E caca o SS = 
(stb)(s2+40s+2500) 2 7 
30.01, 30.1, 30.5, 31 ,35 and 40.0. Which values of b will have minimal effect upon 
the pure second-order transient response? 


4- Given the transfer function G(s) = 


LAB 


1. Using Simulink, add a pole to the second-order system of Prelab la and plot the 
step responses of the system when the higher-order pole is nonexistent, at 
-200, - 24, - 10, and -2. Make your piots on a singie graph, using the Sımuiini 
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LTI Viewer. Normalize all plots to a steady-state value of unity. Record percent 
overshoot, settling time, peak time, and rise time for each response. 


2. Using Simulink, add a zero to the second-order system of Prelab 1 a and plot the 
step responses of the system when the zero is nonexistent, at -200, - 50,-20, - 10, - 5, 
and -2. Make your plots on a single graph, using the Simulink LTI Viewer. 
Normalize all plots to a steady-state value of unity. Record percent overshoot, settling 
time, peak time, and rise time for each response. 


3. Using Simulink and the transfer function of Prelab 3 with a = 3, plot the step 
responses of the system when the value of b is 3, 3.01, 3.1, 3.3, 3.5, and 4.0. Make 


your plots on a single graph using the Simulink LTI Viewer. Record percent 
overshoot, settling time, peak time, and rise time for each 
response. 


4. Using Simulink and the transter function of Prelab 4 with a = 30, plot the step 
responses of the system when the value of b is 30, 30.01, 30.1, 30.5, 31, 35, and 40. 
Make your plots on a single graph, using the Simulink LTI Viewer.Record percent 
overshoot, settling time, peak time, and rise time for each response. 


POSTLAB 


|. Discuss the effeci upon the transient response of the proximity of a higher-order 
pole to the dominant second-order pole pair. 


2. Discuss the effect upon the transient response of the proxim'ty of a zero to the 
dominant second-order pole pair. Explore the relationship between the length of 


the vector from the zero to the dominant pole and the zero's effect upon the pure 
second-order step response. 


3. Discuss the effect of pole-zero cancellation upon the transient response of a 
dominant second-order pole pair. Allude to how close the canceling pole and zero 
should be and the relationships of (1) the distance between them and (2) the 
distance between the zero and the dominant second-order poles. 
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EXPERIMENT [8,4] 
OBJECTIVE 


l- Evaluate the effect of pole and zero location upon the time response of first- and 
second-order systems. 


2- Evaluate the effect of additional poles and zeros upon the time response of 
second-order systems. 


MINIMUM REQUIRED SOFTWARE PACKAGES MATLAB, Simulink, and 
the Control System Toolbox. 


PART 1 
PRELAB 


a : a 3 ; 
l- Given the transfer function G(s) = a evaluate settling time and rise time for 
s+a 
the following values of a: 1, 2, 3, 4. Also, plot the poles. 


2- Giv ransfer function G(s) = =——. - 
2- Given the trans aia 


a- Evaluate percent overshoot, settling time, peak time, and rise time for the 
following values: a = 4, b = 25. Also, plot the poles. 


b- Calcutate the values ot a and b so that the imaginary part of the poles remains 


the same but the real part is increased two times over that of Prelab 2a, and 
repeat Prelab 2a. 


c. Calculate the values of a and b so that the ima 
the same but the 


Prelab 2a. 


ginary part of the poles remains 
real part is decreased (1/2) time over that of Prelab 2a, and repeat 


3. a. For the system of Prelab 2a, calculate the . alues of a and b so that the real part 
of the poles remains the same but the imaginary part is increased two times over that 
of Prelab 2a, and repeat Prelab 2a. 


b. For the system 


0: Prelab 2a, calculate the values of a and b so that the real 
of the poles remains the same but the imaginary part is increased four times 


part 
over that of Prelab 2a. and repeat Prelab 2a. 


Por the svetem af sla}, > Jilat PETE eee ia i ' > 
4. a. For ihe system: of P-elab 2a, calculate the values of a and b so that the damping 


; 
Po 


4: 
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ratio remains the same but the natura! frequency is increased two times over 
that of Prelab 2a, and repeat Prelab 2a. 


b. For the system of Prelab 2a, calculate the values of a and b so that the damping 
ratio remains the same but the natural frequency is increased four times over 
that of Prelab 2a, and repeat Prelab 2a. 


5. Briefly describe the effects on the time response as the poles are changed in each 
of Prelab 2, 3, and 4. 


LAB 


1. Using Simulink, set up the systems of Prelab 1 and plot the step response of 
each of the four transfer functions on a single graph by using the Simulink LTI 


Viewer. Also, record the values of settling time and rise time for each step 
response. 


2. Using Simulink, set up the systems of Prelab 2. Using the Simulink LTI Viewer, 
plot the step response of each of the three transfer functions on 


a single graph. 
Also, record the values of percent overshoot, settling time, peak time, and rise 
time for each step response 


3. Using Simulink, set up the systems of Prelab 2a and Prelab 3. 1 
LTI Viewer, plot the step response of each of 1 
single graph. Also, record the values of perce 
time, and rise time for each step response. 


ing the Simulink 
he three transfer functions on a 
‘vershoot, settling time, peak 


4. Using Simulink, set up the systems of Prelal 2: ind Prelab 4. Using the Simulink 
LTI Viewer, plot the step response of each of the three transfer functions on a 

single graph. Also, record the values of percent overshoot, settling time, peak 

time, and rise time for each step response. 


POSTLAB 


1. For the first-order systems, make a table of calculated and experimental values of 
settling time, rise time, and pole location. 


2. lor the second-order systems of Prelab 2, make a table of calculated and 


experimental values of percent overshoot, settling time, peak time, rise time, and pole 
location. 


3. For the second-order systems of Prelab 2a and Prelab 3, make a table of 
calculated and experimental values of percent overshoot, setiling time, peak 
time, rise ume, and pole location. 
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4. For the s 
calculated 
tim 


=cond-order systems of Prelab 2a and Prelab 4, make a table of 
and experimental values of percent overshoot, settling time, peak 
“> PISE time, and pole location. 


5. Discuss the effects of pole location upon the time response for both first- and 
second-order systems. Discuss any discrepancies between your calculated and 
experimental values. 


PART 2 
PRELAB 


l. a. Given the transfer function G(s) = evaluate the percent over- 


s2+45+25 
shoot, settling time, peak time, and rise time. Also, plot the poles. 


b. Add a pole at -200 to the system of Prelab la. Estimate whether the transient 
response in Prelab la will be appreciably affected. 


c. Repeat Prelab lb with the pole successively placed at -20, - 10, and -2. 


2. A zero is added to the system of Prelab la at -200 and then moved to -50, - 20, 
- 10, - 5, and -2. List the values of zero location in the order of the greatest to the 
least effect upon the pure second-order transient response. 


25b 
(tra) 
(st+b)(s*+4s+25)’ 
3.3, 3.5, and 4.0. Which values of b will have minimal effect upon the pure second- 


order transient response? 


3- Given the transfer function G(s) = let a= 3 andb=3.01, 3.1, 





EDO +a) 

(s+b)(s?+40s+2500) “t ?=30andb= 
30.01, 30.1, 30.5, 31 , 35 and 40.0. Which values of b will have minimal effect upon 
the pure second-order transient response? 





4- Given the transfer function G(s) = 


LAB 


1. Using Simulink, add a pole to the second-order system of Prelab la and plot the 
step responses of the system when the higher-order pole is nonexistent, at 
FIAR 


Ar CD 2 Wats si late Amas : O ey 
200, - 20, - 19, and -2 Make your piots oi a single graph, using the Simulink 
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LTI Viewer. Normalize all plots to a steady-state value of unity. Record percent 
overs':oot, settling time, peak time, and rise time for each response. 


2. Using Simulink, add a zero to the second-order system of Prelab | a and plot the 
step responses of the system when the zero is nonexistent, at -200, a 50,-20, - 10, - 5, 
and -2. Make your plots on a single graph, using the Simulink LTI Viewer. 
Normalize all plots to a steady-state value of unity. Record percent overshoot, settling 
time, peak time, and rise time for each response. 


3. Using Simulink and the transfer function of Prelab 3 with a = 3, plot the step 
responses of the system when the value of b is 3, 3.01, 3.1, 3.3, 3.5, and 4.0. Make 
your plots on a single graph using the Simulink LTI Viewer. Record percent 


overshoot, settling time, peak time, and rise time for each 
response. 


4. Using Simulink and the transfer function of Prelab 4 with a = 30, plot the step 
responses of the system when the value of b is 30, 30.01, 30.1, 30.5, 31, 35, and 40. 
Make your plots on a single graph, using the Simulink LTI Viewer.Record percent 
overshoot, settling time, peak time, and rise time for each response. 


POSTLAB 


1. Discuss the effect upon the transient response of the proximity of a higher-order 
pole to the dominant second-order pole pair. 


2. Discuss th» effect upon the transient response of the proximity of a zero to the 
dominant second-order pole pair. Explore the relationship between the length of 


the vector from the zero to the dominant pole and the zero's effect upon the pure 
second-order step response. 


3. Discuss the effect of pole-zero cancellation upon the transient response of a 
dominant second-order pole pair. Allude to how close the canceling pole and zero 
should be and the relationships of (1) the distance between them and (2) the 
distance between the zero and the dominant second-order poles. 
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